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ABSTRACT 

We investigate the effects of a variety of ingredients that must enter into a realistic model 
for disk-galaxy formation, focusing primarily on the Tully-Fisher (TF) relation and its scat- 
ter in several wavebands. In particular, we employ analytic distributions for halo-formation 
redshifts and halo spins, empirical star-formation rates and initial mass functions, realistic 
stellar populations, and chemical evolution of the gas. Our main findings are: (a) the slope, 
normalization, and scatter of the TF relation across various wavebands is determined both by 
halo properties and star formation in the disk; (b) TF scatter owes primarily to the spread in 
formation redshifts. The scatter can be measurably reduced by chemical evolution, and also in 
some cases by the weak anti-correlation between peak height and spin; (c) multi-wavelength 
constraints can be important in distinguishing between models which appear to fit the TF 
relation in / or K; (d) successful models seem to require that the bulk of disk formation can- 
not occur too early (z > 2) or too late (z < 0.5), and are inconsistent with high values of 
Flo; (e) a realistic model with the above ingredients can reasonably reproduce the observed 
z = TF relation in all bands (B, R, I, and K). It can also account for the z — 1 _B-band 
TF relation and yield rough agreement with the local B and K luminosity functions and B- 
band surface-brightness-magnitude relation. In such a model, the near-infrared TF relation at 
z = 1 is similar to that at z — 0, while bluer bands show a markedly steeper TF slope at high 
redshift. The remarkable agreement with observations suggests that the amount of gas that is 
expelled or poured into a disk galaxy must be small (though small fluctuations might serve to 
better align _B-band predictions with observations), and that the specific angular momentum 
of the baryons must roughly equal that of the halo; there is little room for angular momentum 
transfer. In an appendix we present analytic fits to stellar-population synthesis models. 

Key words: cosmology: theory — galaxies: formation — galaxies: spiral — galaxies: kinemat- 
ics and dynamics. 



1 INTRODUCTION 

Spiral galaxies are particularly important in the study of galaxy for- 
mation, as they are believed to undergo a relatively smooth for- 
mation process, and serve as the building blocks in the formation 
of other galactic systems through mergers. Thus, spiral galaxies 
should be the easiest to model and should provide clues as to the ba- 
sic physics underlying galaxy formation. Various lines of observa- 
tional evidence serve to guide our understanding of spirals, includ- 
ing their measured luminosity function (LF), surface-brightness 
distribution, star-formation history, chemical composition, and dy- 
namical properties. Of particular significance is the Tully-Fisher 
(TF) relation, a remarkably tight correlation between the luminos- 
ity and rotation speed of spirals. For a given wavelength, A, the TF 
relation obeys the form L\ = A\ V c lx , or 

M A = a A + &AlogV;, (1) 



where M\ is the absolute magnitudejj, and 6a = —2.57a is the 
slope of the relation. The debate continues as to whether this re- 
lationship results primarily from initial conditions, i.e., from the 
properties of the parent halo (Dalcanton, Spergel, & Summers 
1997; Mo, Mao, & White 1998; Firmani & Avila-Reese 1998; 
Avila-Reese, Firmani, & Hernandez 1998; Navarro & Steinmetz 
2000; Mo & Mao 2000), self-regulating feedback processes associ- 
ated with star formation in the disk (Silk 1997), or a combination of 
both (Heavens & Jimenez 1999; Somerville & Primack 1999; van 
den Bosch 2000). 

A TF relation arises quite naturally if one simply assumes that 
galactic halos formed at roughly the same time and that luminosity 



* A subscripted M is understood to refer to absolute magnitude in a given 
band, while an M without a subscript refers to mass. 
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is proportional to the baryonic mass, which is in turn proportional 
to the halo mass. More realistically, the luminosity may depend on 
the galactic spin, as disks formed in high-spin halos will be larger 
and thinner, and thus lead to lower star-formation rates. Still, even 
with spin, a TF relation arises if halos all formed at roughly the 
same epoch. In both cases (with and without spin), scatter in the 
redshift of halo formation should lead to scatter in the central den- 
sities of the halos and thus to scatter in the TF relation. 

Several groups have made considerable progress in under- 
standing spiral-galaxy properties along these lines using semi- 
analytical models (SAMs). Eisenstein & Loeb (1996, hereafter 
EL96) used Monte-Carlo realizations of halo-formation histories 
to calculate the minimum TF scatter which should arise from the 
spread in halo-formation times. They concluded that unless spirals 
form at z <; 1, without subsequently accreting much mass, the TF 
relation cannot arise simply from initial conditions, but must in- 
stead owe to some feedback mechanism which decouples the lumi- 
nosity from the halo history. 

More recently, other groups have investigated detailed SAMs 
of disk-galaxy formation which incorporate such features as forma- 
tion histories derived directly from A^-body simulations, universal 
halo profiles, adiabatic disk contraction, bulge formation via stabil- 
ity criteria, star formation, supernova feedback, dust, cooling, and 
mergers (Somerville & Primack 1999; van den Bosch 2000; Fir- 
mani & Avila-Reese 1998; Navarro & Steinmetz 2000). Their con- 
clusions differ as to the relative importance of initial conditions vs. 
feedback from star formation and/or supernovae in defining the TF 
relation. These studies generally agree, however, that the spread in 
halo-formation redshifts is a significant source of scatter in the TF 
relation and that reconciling models with TF observations seems to 
require a low matter density (fio ~ 0.3 - 0.5) and disk formation 
at high redshift (z <; 1). 

One powerful test of such models, which has not been fully ap- 
preciated, is the simultaneous comparison of their predictions to TF 
data from several wavebands. Some previous work has considered 
only an assumed value for the mass-to-light ratio in a given band, 
rather than using stellar population models to predict broad-band 
magnitudes, and many authors have investigated TF predictions for 
only a single waveband, typically / or K, where the observed TF 
scatter is the smallest [~ 0.4 mag in the most carefully defined 
samples (Willick et al. 1995, 1996, 1997; Tully et al. 1998)]. Since 
these wavelengths measure primarily the oldest, shell-burning, stel- 
lar populations, such measurements are, by construction, effec- 
tively probing only the total mass of the galaxy. Thus, while some 
authors have claimed success in fitting the near-IR TF relation, their 
predictions in bluer bands, where many model ingredients would be 
most strongly manifested owing to the younger populations probed, 
have gone unchecked. With realistic stellar-population models, the 
leverage gained by spanning a range of wavelengths should there- 
fore be crucial in distinguishing models which produce similar 
near-infrared TF predictions, and help to assess the importance 
of various SAM features. Heavens and Jimenez (1999; hereafter 
HJ99), used a simple halo and disk model combined with empiri- 
cal star-formation properties to investigate the role of star forma- 
tion in the TF relation, did examine the simultaneous constraints 
from various wavebands, but did not address some of the other key 
considerations listed above. 

Our approach here will focus on constructing, as much as pos- 
sible, a model "from the ground up," i.e., starting with a minimal 
set of simple, well-motivated assumptions, and individually inves- 
tigating the impact of various plausible modifications. In essence, 
we are asking whether a plausible model with fewer free parame- 



ters^ than are usually incorporated into current SAMs can reason- 
ably pass various observational tests. Though the resulting model 
will lack many of the sophisticated features of other SAMs and will 
surely be an oversimplification in many respects, it may shed im- 
portant light on the issue of which physical ingredients are truly 
essential in determining certain disk properties. 

We build on the work of HJ99, using their halo/disk model 
and stellar populations code, but generalizing it to include a vari- 
ety of important features, such as analytic models for the distribu- 
tion of halo-formation redshifts and spins (as well as the predicted 
anti-correlation between the two). This allows us to assess not only 
the impact of these distributions on TF predictions, but the depen- 
dence on cosmological parameters and the power spectrum as well. 
We also include for the first time an analytic description of chem- 
ical evolution. This measurably reduces the scatter in the TF rela- 
tion, which is found to arise mainly from the spread in formation 
times. The overall result will be a model for disk galaxies which 
passes a host of observational tests, reproducing the observed TF 
relation, locally and at z = 1, in all relevant wavebands with the 
right magnitude of scatter. It also succeeds in roughly fitting the 
observed surface-brightness distribution of spirals and, within its 
limitations, reasonably reproduces the observed B- and AT-band 
luminosity functions. The model necessitates cosmological param- 
eter values which are in line with current estimates, and requires 
that the bulk of disk formation occur in the range 0.5 < z < 2. 
This agreement implies that, although plausible, other mechanisms 
which would add or remove gas from the disk (such as mergers, gas 
expulsion by supernovae, etc.), and thereby distort the predicted lu- 
minosity, are not necessarily required to successfully predict these 
global disk-galaxy properties. 



2 THE MODEL 

2.1 Our Starting Point 

We build on the skeletal disk-galaxy model of HJ99. We review 
the basic ingredients here and refer the reader to HJ99 for a more 
complete description. This model assumes that after dark halos 
separate from the Hubble flow and collapse, they subsequently 
relax to an isothermal sphere of mass M, and the baryons in- 
stantaneously settle into a central disk with mass Md = rridM, 
where m c i = Q^/Qo, Qo is the present-day non-relativistic mat- 
ter density in units of the critical density, and the baryon den- 
sity, Sib, is set by nucleosynthesis, Qb = 0.019fo -2 (Tytler et 
al. 1999). The disk is assumed to have an exponential profile, 
E(i?) = Eo exp(— R/Rd), with central surface density Eo and 
scale length Rd, such that Md = 27rEoi?d 2 . Since the assumptions 
of a singular isothermal sphere and instantaneous disk formation 
cannot strictly be true, the apparent success of the model will argue 
against any strong influence of details such as the profile shape (Mo 
& Mao 2000) or disk-contraction process on the resultant proper- 
ties of spirals. 

With the gaseous disk in place at the time of formation, star 
formation proceeds according to an empirical Schmidt law (Ken- 
nicutt 1998), such that the star-formation rate (SFR) depends only 
on the local gas surface density. Given the conservation of baryonic 

t Most of the freedom in our models comes from the cosmological- 
parameter choices, which we generally fix to be consistent with current es- 
timates; other choices for these parameters would yield other TF relations, 
usually inconsistent with the data. 
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mass, the total SFR, as well as remaining gas fraction, can then be 
obtained as functions of time since formation by integrating over 
the disk [see equations (3)-(9) in HJ99]. The model ignores any 
gas returned to the ISM by stars as well as late infall of fresh gas. 

Using the SFR and an assumed stellar initial mass function 
(IMF), the properties of the stellar populations formed at each 
time interval are calculated using the spectrophotometric stellar- 
evolution code of Jimenez et al. (1998, 2000), and broad-band 
magnitudes derived from the resulting spectra are integrated over 
the life of the galaxy to yield luminosities in various bands. Other 
stellar-population codes generally yield results which agree to 
within 0.1-0.2 mag and have comparable dispersions (Jimenez et 
al. 2000). Since we are integrating the SFR over time, weighted 
by the specific luminosities of model stellar populations as a func- 
tion of age, metallicity, spin, etc., the model can predict disk lu- 
minosities and surface brightnesses (in a variety of wavebands) at 
any point in their evolution. More details about the calculation of 
the luminosities of the stellar populations, as well as some possibly 
useful analytic fitting formulas, are provided in an Appendix. 

To summarize, for a given set of cosmological parameters, the 
model of HJ99 takes as input a total galactic mass, M, and spin, 
A, and can output, as functions of time, the disk luminosity and 
surface brightness in various wavebands, as well as the metallicity. 

2.2 Spherical-Collapse Model 

We will be interested in the relation between the size and mass of a 
galactic halo, and this is fixed by the formation redshift — halos of a 
given mass that undergo collapse earlier are expected to be smaller 
and denser — through the spherical-collapse model. HJ99 used a 
simple approximation for the relation between the mass, size, and 
collapse redshift, neglecting the full impact of a possible cosmolog- 
ical constant. Since our aim is to carefully study the effects of vari- 
ations in the formation-redshift distribution, we treat the spherical- 
collapse model more precisely, following the approach of Wang & 
Steinhardt(1998). 

The model will of course depend on the background cos- 
mogony, which effectively sets the initial conditions and dictates 
the course of structure formation. We consider various cold-dark- 
matter (CDM) models with a power spectrum given by P(k, z) = 
AD 2 (z)k n T 2 (k), where A is the overall amplitude, D(z) is the 
linear growth factor, and k n is the primordial power spectrum (we 
shall assume an untilted, n = 1 primordial spectrum throughout). 
For the transfer function, T(k), we adopt the functional form of 
Bardeen et al. (1986), parameterized by the shape parameter T. 

We restrict out study to flat geometries (Qo + £7a = 1, where 
SIa is the energy density associated with the cosmological con- 
stant, in units of the critical density), and consider various values 
of h, the Hubble parameter (Ho) in units of 100 km s _1 Mpc -1 . 
The presence of a cosmological constant introduces several mod- 
ifications to the standard spherical-collapse model derived for the 
Einstein-de Sitter case, which are often overlooked. These include 
a nonzero surface-pressure term in the virial relation, which in turn 
alters the mass-circular velocity relationship, and introduces fio- 
dependent corrections to p co u — the density of virialized objects, 
in units of the present background density, given by a value of 
~ 179 at z = in the Einstein-de Sitter (EdS) case — and to 5 C , 
the overdensity (S = p/p — 1) threshold for collapse, typically 
taken to be the linearly-extrapolated value of ~ 1.69 in the EdS 
case (see equations (2), (5), (6), and (7) and related discussion in 
Wang & Steinhardt 1998). We do not reproduce all of these results, 
but merely write the equation for the quantity which interests us di- 



rectly, namely the circular velocity of the disk, which can be related 
to the virial velocity of particles in the halo: 



Vc 



2 -— 

3 v vlI 



(2) 



GM# \M>A c /2 

2 1 - n 



2/3 



(l + z) X 



where A c = p co u / pb is the fio-dependent ratio of the halo to back- 
ground density at collapse. We stress that these results can, in gen- 
eral, differ substantially from their EdS counterparts. 

We incorporate these ACDM generalizations into the Press- 
Schechter (PS) theory, which will be our starting point for calcu- 
lating halo abundances and formation times. The theory states the 
familiar result that the comoving number density of halos in a mass 
interval dM about M at redshift z is given by 



n(M,z)dM = - 



2 p b 



da(R,z) 



ttM a 2 (R,z) 



dM 



(3) 



exp 



2a 2 (R,z) 



dM, 



where i? 3 = 3M/4ivpb, pb is the constant, comoving background 
matter density, o(R, z) is the rms fractional density perturbation in 
spheres of radius R, and 5 C is here understood to depend on redshift 
via an empirically derived relation (L. Wang, private communica- 
tion). 



2.3 Formation-Redshift Distribution 

In standard hierarchical models of structure formation, halos of a 
fixed mass form over a range of redshifts. HJ99 examined only 
fixed values for the redshift of formation, z/, for halos, but con- 
sidered several values to illustrate the effect of a spread in z/. 
EL96 performed Monte-Carlo realizations of halo formation his- 
tories, using the merger- tree approach (Lacey & Cole 1993, 1994) 
and a spherical-accretion model, and concluded that the minimum 
TF scatter resulting solely from the calculated spread in halo for- 
mation redshifts is already uncomfortably larger than that observed. 
They concluded that satisfying the upper limit of a ~ 10% relative 
error in velocity dispersion requires that Az//(1 + z/) <, 20%. 
Other workers (e.g., Bullock et al. 1999) have looked at the scatter 
in formation redshifts (or the closely-related scatter in "concentra- 
tion parameters") obtained from numerical simulations and found 
better agreement with the observed scatter. 

Although numerical simulations should in principle provide 
the best way to evaluate the formation-redshift distribution, these 
have limitations in practice. In particular, it is difficult to collect 
enough statistics to determine the distributions for different masses 
and/or to determine how these distributions depend on the cosmo- 
logical parameters or the power spectrum. As the analytic distri- 
butions discussed below demonstrate, it is indeed to be expected 
that the formation-redshift distribution (and thus resulting scatter in 
the TF relation) should depend on the mass, cosmological param- 
eters, and power spectrum. Thus, although they can provide some 
order-of-magnitude estimates, results from numerical simulations 
can be mis-applied in SAMs if they are determined, say, for one 
mass from a simulation with a particular choice of cosmogony, and 
then applied to other masses and/or cosmogonies. Moreover, ana- 
lytic methods allow for more obvious and direct insight into the de- 
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pendence of the model on various parameter choices and assump- 
tions. 

Therefore, we follow Viana & Liddle (1996), and explore two 
plausible analytic models for the distribution of halo formation red- 
shifts. The first, denoted as the 'S' distribution, is that of Sasaki 
(1994), which simply uses the PS formalism to calculate the for- 
mation rate of bound objects, weighted by the probability of arriv- 
ing at some later time without merging, under the assumption that 
the destruction efficiency is independent of mass. In this model, the 
distribution of formation redshifts of halos with mass M is given 
by 



dns 8 C 2 n(M, Zf) da(R, Zf) 

dzf a 2 (R,zf) a(R,0) dzf 



(4) 



The S distribution has the advantage of being independent of how 
one defines a new halo, but the assumption of self-similar merging 
might be questioned (see, e.g., Percival, Miller & Peacock 2000). 

The second distribution, denoted as the 'LC distribution, ad- 
dressed this shortcoming by employing the merger-tree approach 
of Lacey & Cole (1993, 1994). In this formalism, halos continu- 
ally grow with time, and so one must explicitly define the time at 
which a particular halo has come into existence^] (e.g., the epoch of 
mass doubling). Lacey & Cole (1993) derived expressions for the 
expected z / distribution of halos both using analytic counting ar- 
guments and a Monte-Carlo approach to generate merger histories 
for M(t) explicitly. They found that the former approach yielded 
excellent agreement with iV-body results, while the latter overesti- 
mated halo ages (Lacey & Cole 1994). Their results did not depend 
strongly on the value of the local slope of the power spectrum. The 
z f distribution in this model is given by 



driLC 
dzt 



where 



-p(w(z f )) 



dw(zf) 
dzt 



(5) 



p(w(zf)) 



w(z f ) = 



2w{zf){f- - l)erfc 



V2 



-IT 1 -2) 



exp 



5c(a(M,a)/a(M,z f )-l) 

y/<7*(fM,0)-<T*(M,0) 



(6) 



(7) 



and / is the fraction of halo mass assembled by formation red- 
shift Zf. Thus at fixed redshift, models with higher values of / 
produce younger galaxies. Lacey & Cole adopted / = 0.5, iden- 
tifying the mass doubling epoch as the nominal point at which ha- 
los are born, while Viana & Liddle (1996), looking specifically at 
the mass-temperature relation for galaxy clusters, use a value of 
/ = 0.75. Note that in any of these models, the dependence on the 
power spectrum is only through a(M), the rms fractional density 
perturbation for the mass scale M; the dependence on fio and h 



t This ambiguity arises because in the excursion-set formalism, one fol- 
lows a random walk of trajectories of the spatially-filtered <5 to the first 
up-crossing above some threshold value. The "tagged" mass is statistically 
defined as the expectation value for the mass of a halo in which a given 
particle will end up, but this will in general differ from the actual mass in a 
given realization. This ambiguity can lead to predictions of negative proba- 
bility densities at z ~ in models with primordial power-spectrum indices 
n > 0. 



i — 1 — 1 — 1 — 1 — i — 1 — 1 — 1 — 

S 

LC, f=0.5 

LC, f=0.75 

- LC, f=0.9 




Figure 1. Formation-redshift distribution for the S distribution (solid line) 
and for the LC distribution with / = 0.5, 0.75, and 0.9 (dotted, short- 
dashed, and long-dashed lines, respectively), assuming a ACDM model 
with Q = 0.3 and h = 0.65. Note the appreciable width of the S dis- 
tribution, and the decreasing width and peak z e values for LC distributions 
with increasing /. 



enter through the linear-theory growth factor and in the spherical- 
collapse physics, and through the star-formation model by fixing 
ages, baryon fraction, etc. 

Figure [l] depicts the distribution of formation redshifts for ha- 
los of 1O 12 M0 both for the S distribution, and for the LC distri- 
bution with / values of 0.5, 0.75, and 0.9 for comparison. These 
results assume a COBE-normalized ACDM model with fi = 0.3 
and h — 0.65 (yielding, in this case, a mean peak height of 
v = 0.69, and a s = 1.08). Note the large width of the S distribu- 
tion, peaking at z ~ 1.0, but with a significant tail extending out to 
z ~ 4—5, and appreciable ongoing formation today^} By contrast, 
LC distributions with larger values of / produce much narrower 
distributions which peak at lower redshifts and fall more rapidly to 
high z. 



2.4 Spin Distributions 

Assuming constant specific angular momentum, the disk scale 
length can be related to A, the spin parameter of the halo (e.g., Mo, 
Mao, & White 1998), 



Rd 



XGM 



(8) 



fixing the initial surface density of the gas and thus affecting the 
SFR, luminosity, and surface brightness. Therefore, it might be sup- 
posed that the predicted TF relation and other properties would de- 
pend sensitively on spin. HJ99 considered only fixed values for A, 
rather than accounting for its detailed distribution. The spin distri- 
bution is usually taken to be a log-normal centered at A ~ 0.5, 



S This long tail will disappear if one uses the formalism developed in Per- 
cival, Miller & Peacock 2000 
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as indicated by numerical simulations. However, as is the case for 
the formation-redshift distribution, heuristic arguments suggest that 
the distribution of spin parameters should depend on the mass, the 
formation redshift, and possibly on the cosmological parameters 
and/or power spectrum. If so, the effect of scatter in the spin pa- 
rameter on the TF relation could be different than that inferred from 
numerical simulations. 

Heavens & Peacock (1988) calculated the distribution of tidal 
torques acting on matter in the vicinity of mildly nonlinear density 
maxima, assuming a spherical-accretion model to calculate bind- 
ing energies, and derived the resulting spin-parameter distribution. 
They found an anti-correlation between peak height and spin pa- 
rameter, but pointed out that the intrinsically broad range of the A 
distribution swamped the systematic shift with peak height, v, re- 
sulting in a fairly weak anti-correlation. Catelan & Theuns (1996) 
extended these results to obtain the joint PDF for spin parameters 
and peak masses using the distribution of peak shapes in different 
CDM models. They confirmed the broadness of the spin-parameter 
distribution, and the anti-correlation with peak height (which owes 
essentially to the fact that higher peaks will generally be more 
spherical, and thus harder to spin up). Their results, however, ex- 
hibit a systematic shift toward higher A values than those of Heav- 
ens & Peacock, such that v = 1 peaks have rms A values of ~ 0.15, 
instead of 0.05, thus closer to the observed value of 0.5 for spirals. 
The anticorrelation conjecture has not been fully tested by simu- 
lations, but some support comes from Ueda et al. (1994). Lemson 
& Kauffmann (1999) argued against an environment-dependent A. 
However, the environment was defined on a 10 h^ 1 Mpc scale, 
and we do not expect strong correlation between the density field 
smoothed on galaxy scales and this large scale. 

Since the indications for a distribution peaked at higher spins 
and an anti-correlation with peak height are good, but not yet well- 
established, we consider two models, with and without the joint 
distribution in A and v. For the former, we adopt the joint distribu- 
tion function (Catelan & Theuns 1996), 



P{\\v) = ^ 



0.68 
A 



1 + 0.02 



A 



Ao(k) 



exp 



log 2 (A/Ao(Q) 
0.98 



(9) 



where Ao(f) = O.lli^ 1 ' 1 and v = u(z) = 8 c (z)/a(R, z). Since, 
at a fixed epoch, there is not a one-to-one relationship between v 
(defined for the spatially averaged overdensity field) and mass, ow- 
ing to the distribution of halo shapes, equation (12) effectively cor- 
responds to an average of v over all shapes of a given mass.p| 

2.5 Chemical Evolution 

HJ99 assumed the stellar populations to have constant metallicity 
fixed at the solar value for all time. More realistic models should, 
of course, account for chemical evolution. To study the effect of 

^ We point out that the binding energy used to calculate A is derived using 
a spherical-accretion model, but that for the lower, more irregular peaks, 
this is likely only a lower limit to the actual binding energy, resulting in an 
underestimate of A. Combined with the arguments above from Catelan & 
Theuns (1996), one can push the predicted spin values for spirals (if these 
are indeed associated with lower-y peaks, as opposed to ellipticals, which 
are associated with higher peaks) closer to the observed values, such that 
one needn't rely as heavily on dissipation to spin up spirals. 



chemical evolution on the TF relation, we assume that galaxies are 
(chemically) closed boxes, for which analytic results for the evolu- 
tion of metallicity exist (e.g., Pagel 1997). This assumption seems 
to be justified by the detailed hydrodynamical computations of a 
multi-phase ISM by Mac Low & Ferrara (1999). In order to re- 
lax the instantaneous recycling approximation we adopt a value of 
0.03 for the yield in order to account for the delayed production 
of Fe by Type-la supernovae. We also set the zero-point redshift 
dependence so as to comply with observations of damped Lyman- 
a systems, which suggest that the ISM is already enriched up to 
a certain metallicity at a given redshift. Adopting the relation de- 
rived by Pettini et al. (1999; see Figure 8 therein), we have that the 
metallicity, Z, evolves with redshift, z, as 



Z{z,z f ) = 0.03Z© In 



0.28Z Q 10"°- 25z , (10) 



where the solar metallicity Zq = 0.02 and f g (z, Zf) is gas fraction 
of the disk. It will be shown that chemical evolution will be an 
important factor in reducing the predicted scatter in the TF relation. 



3 RESULTS 

We first study the respective impacts of the various ingredients 
which have entered into our galaxy-formation model. In Figures 
^| to [l2|, we illustrate examples of how our predictions for the TF 
relation in the B, R, I, and K bands change as the model inputs 
are varied. These results will allow the reader to infer how the pre- 
dictions of any of our models, or any other models that appear in 
the literature, would change with different input physics or assump- 
tions. In each graph, the solid line delineates the power-law fit to the 
TF-relation prediction (along with 1-cr errors given by the dashed 
lines, and denoted in each plot by 'cr'), and the four scattered-dot 
curves trace the spread in the predicted TF relation for four fixed 
masses (10 10 , 10 11 , 10 12 , and 1O 13 M ), using ~ 150 points each. 

Overplotted on each graph are open symbols representing data 
from Tully et al. (1998). These data, comprised of spiral galaxies 
in the loose clusters of Ursa Major and Pisces, represent one of 
the most carefully defined TF samples. The inclination-corrected 
FWHM of the lines has been converted to W l R , which approxi- 
mates to 2V"c, and the data have been corrected for dust extinction. 
The scatter in the B, R, I, and K bands is 0.50, 0.41, 0.40, and 
0.41 mag, respectively. While many authors compare their models 
to TF data by defining a "fundamental" TF relation derived from 
the particular data set in question, we argue that, given the obser- 
vational selection effects inherent in defining such samples, there 
is no particular significance or accuracy associated with any one 
determination of the slope and normalization of the TF relation. In- 
stead, we examine how well these data fit to our model, in a \ 2 
sense, using the measured scatter as an estimate of the error associ- 
ated with each point. For each plot we calculate the resulting \ 2 and 
denote by 'p' the probability of obtaining this large a value of x 2 
given that the model is correct; values of p below 10~ 3 are listed as 
zero. In addition, since the data have excluded spirals which show 
evidence of merger activity or disruption in the form of starbursts, 
we exclude from our model galaxies those having B — R < 0.3. 



3.1 Formation-Redshift Distribution 
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3.1.1 S Distribution 

As discussed above, the distribution of formation redshifts should 
lead to some scatter about the TF relation, and Figure ^| illustrates 
this effect using the S distribution. We show results for two cos- 
mogonies to indicate the dependence of the results on some of the 
cosmological parameters. The left-hand panels of Figure ^ depict 
the results for the S distribution of zj in an EdS universe, while the 
right-hand panels show the results for a fiat, A-dominated universe. 
Here, and unless otherwise indicated, we shall take EdS models 
to have h = 0.65, a power-spectrum shape parameter^ V = 0.2, 
and an rms density contrast fluctuation over 8/i _1 Mpc spheres of 
ag = 0.5. For ACDM models, we assume Qq = 0.3, h = 0.65, 
a COBE-normalized power spectrum (Bunn & White 1997) with 
T = floh = 0.195, and for all cosmogonies we assume A = 0.05, 
a fixed metallicity given by the solar value (Z — Zq), and a 
Salpeter IMF, d log N* /d log M* = -a with a = 1.35, all un- 
less otherwise indicated. 

To understand the results, we begin by examining the predic- 
tions for the if-band TF relation in the ACDM model. Here we find 
that the spread in Zf translates directly into a spread in V c (with ear- 
lier formation implying higher circular velocity) which widens for 
lower masses, since these can form over a broader range of red- 
shifts stretching back to higher z. This broadening of the TF with 
lower V c arises generically from a realistic formation-distribution, 
and it is observed in most TF samples. The predicted Mk, how- 
ever, varies very little for galaxies of a given mass, since this wave- 
length primarily traces the older stellar population. In general for 
our models, the distribution in Zf is found to be the primary source 
of scatter in the TF relation. 

As we go to bluer bands, we find that for a given mass, 
younger galaxies become brighter, increasing the spread in the TF 
relation. This is a simple consequence of equation (8) in HJ99, 
which states that the SFR is a rapidly declining function of age. 
There is a competing effect, however, from the fact that galaxies 
that are too young will not have had enough time to convert much 
gas into stars, and so will become dimmer. For this reason, our 
plots often outline a peak magnitude for galaxies of a given mass, 
in those cases where the Zf distribution produces a large fraction of 
very young galaxies. Such galaxies would likely exhibit the prop- 
erties of starbursts and be excluded from carefully-defined TF sam- 
ples. Indeed, model galaxies removed by our color selection crite- 
rion (B — R < 0.3) tend to be those on the younger side of these 
peaks. 

Overall, these S distribution predictions do not fit the TF data 
in any waveband, as the models are too faint in every case, and the 
predicted scatter here is too large by about 0.7 mag in B and 0.3 
mag in K. The problem owes largely to the fact that the S distri- 
bution has a substantial high-redshift tail producing older, fainter 
galaxies with high V c . Thus, if we assume star formation proceeds 
according to a Schmidt law, TF observations imply that the ma- 
jority of present-day disks could not have formed much beyond a 
redshift of 2-3. 



II Since values of Qoh ^> 0.3 already seem to be ruled out (e.g., Peacock 
& Dodds 1994), V in the case of EdS models is defined simply as a fitting 
parameter in the CDM transfer function. In the case of ACDM models, we 
formally take T = QqH. 



3.1.2 LC Distribution 

Figure ^ depicts corresponding results for the same parameter 
choices, but using the LC distribution with / = 0.75. This distri- 
bution has both a lower peak value of Zf and a smaller FWHM than 
the S distribution (see Figure |l|), producing significantly less scat- 
ter in the predicted TF relations and a larger fraction of younger 
galaxies, as compared with Figure ^. Despite the smaller scatter, 
these models still do not fit the observed TF relation in the vari- 
ous bands. Note that the EdS model here actually yields excellent 
agreement in B, but predicts galaxies which are too faint in I and 
K. 

In Figures ^ and ^ we examine LC distributions with values 
of / = 0.5 and / = 0.9, respectively. Lacey & Cole (1993, 1994) 
found that the statistics of the time at which present-day halos ac- 
creted half their mass (i.e., / = 0.5) provides an excellent fit to 
the results of iV-body simulations. Here we see that the / = 0.5 
ACDM model provides an excellent fit to the K-, R-, and /-band 
TF relations but predicts galaxies which are too faint in B, thus 
highlighting the importance of using multi-wavelength constraints. 
This model also produces TF scatter which is about 0. 1 to 0.4 mag 
too large, going from K to B. The EdS predictions here have less 
scatter, but are again too faint in every case. In general, we find 
that LC distributions with plausible values of / fare better than S 
distributions at matching TF data, suggesting disk formation in the 
range < z < 2 (see below). 

The / = 0.9 models have a Zf distribution which produces 
galaxies only in a narrow interval centered at low redshift (see Fig- 
ure [l]). This results in extremely small TF scatter in part because ef- 
fectively no galaxies form above a relatively low maximum redshift 
of z ~ 1, but also because many of these recently formed galax- 
ies fall on the young side of the magnitude peak and are excluded 
by our color selection. In effect, this zj distribution has "pushed" 
galaxies over to the faint side of the magnitude peak, where their 
colors resemble those of starbursts, and the low scatter owes to the 
small number of galaxies "older" than the peak, particularly in K. 
Thus, while the intrinsic TF scatter of all galaxies in this model is 
high, the observed scatter, based on observational selection, can be 
quite low, and this effect is important when considering recently 
formed populations. 

Irrespective of the scatter, however, the / = 0.9 models 
fail to fit the TF slope and normalization. For the reasonable 
cosmological-parameter values in the right panels, the predicted 
galaxies are very young and thus excessively bright. Thus, TF con- 
straints also imply that the majority of disk formation cannot have 
happened at z < 1. Of course, from a physical standpoint, such 
large values of / are anyway unreasonable. Note that even with 
/ = 0.9, the EdS model still produces galaxies which are too faint 
in K, despite yielding excellent agreement in /. 

3.2 Dependence on Cosmogony 

3.2.1 Power-Spectrum Shape and Amplitude 

In Figure |^ we explore the impact of individually changing the 
shape and normalization of the power spectrum, using the S dis- 
tribution in an EdS universe. In the left (right) panels we adopt 
T = 0.5 (T = 0.2) and erg = 0.5 (erg = 1.0). In both cases, 
the resulting galaxies are several magnitudes too faint and produce 
dramatically larger scatter in their TF relation. In the case of higher 
as, halos of a given mass will be smaller, and correspond to rarer 
peaks. These collapse at earlier times, leading to larger V c and 
older, dimmer galaxies. Extending the range of collapse redshifts 
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Figure 2. S-distribution predictions. The left panels show the predictions for the present-day TF relation in B, R, I, and K for the S distribution assuming an 
EdS universe with h = 0.65 and a CDM power spectrum with T = 0.2 and erg = 0.5. The right panels display the results for a COBE-normalized ACDM 
model with Qo = 0.3, h = 0.65, and T = Qoh. Both assume a Salpeter IMF, a constant spin parameter of A = 0.05, and constant metallicity fixed at the 
solar value. The model predictions for halos of mass 10 10 , 10 11 , 10 12 , and 10 13 Mq are indicated by the four scatter-point curves in each panel, and the solid 
lines are the fits to these predictions, with 1-ct errors indicated by the dashed lines. The open symbols are data from Tully et al. (1998). Note the EdS model, 
in particular, underpredicts the luminosity in all bands. 
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Figure 3. LC-distribution predictions. Same as Figure El but using the LC distribution with / = 0.75. 
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Figure 4. Same as Figure hi but using / = 0.5. Note how the A model yields excellent agreement in K, R, and I, but provides a relatively poorer fit in 
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Figure 5. Same as Figure but using / = 0.9. In this case, all disks form very recently and many fall on the younger side of the magnitude peak. These 
effects combine to yield very little TF scatter for the few relatively older galaxies which fall red- ward of our color selection criterion. The overall agreement 
however, is poor. 
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Figure 6. Dependence on the power spectrum shape and amplitude. TF predictions for the S distribution in an EdS universe for V = 0.5 and erg = 0.5 (left 
panels), and for V = 0.2 and erg = 1.0 (right panels). Comparing with Figure ^ we find that high T and high erg both lead to earlier collapse, producing 
fainter disks with a much larger TF scatter. 
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Figure 7. Dependence on the Hubble constant. Same as Figure^, but using h = 0.5. The resulting changes in baryon fraction and age lead to a brightening 
of the TF relation, particularly for low Qo- 
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Figure 8. Same as Figure El but using h = 0.75. 
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Figure 9. Effect of joint distribution in spin and peak height. Same as Figure |2| but now incorporating both the predicted shape of the spin-parameter 
distribution and the anti-correlation between spin and peak height. These results differ little from those of Figure [jj the joint distribution in A and v scatters 
the predictions, by a small amount, along the TF relation itself, tightening the TF relation in some cases. 
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Figure 10. Same as Figure W, but now incorporating the joint PDF in A and v into the LC results. 
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Figure 11. Effect of chemical evolution. TF predictions for an Co = 1 S distribution (left panels) and an / = 0.75 LC distribution with Ho = 0.3 (right 
panels), both with evolving metallicity. Note the smaller TF scatter as compared to the corresponding results in Figures Pj and H. 
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Figure 12. Dependence on the IMF. TF predictions for an / = 0.75 LC distribution with Co = 0.3, assuming a Salpeter IMF with index a = 0.95 (left 
panels), and a Scalo IMF (right panels). These IMFs produce relatively more massive stars than the a = 1.35 Salpeter IMF investigated previously, and thus 
brighter TF predictions than those in the right panels of Figure H. 
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back to higher z also substantially increases the scatter. Similarly, 
in these CDM models, increasing Y shifts more power to smaller 
scales (R < 10 Mpc), leading to the same result. Though we do not 
explicitly explore variations with n, it can be inferred that values of 
n > 1, which also tilt power towards smaller scales, will have a 
similar effect. Conversely, values below unity might lead to better 
agreement with the data, but there are already strong constraints on 
the tilt (e.g., Kamionkowski & Buchalter 2000). 



3.2.2 Hubble Constant and Matter Density 

Given the uncertainties in our present understanding of galaxy for- 
mation, it makes little sense to attempt to constrain cosmological 
parameters using TF data. Rather, with these parameter values dic- 
tated by other, more direct tests, we should employ the TF relation 
to gain insight into these details. Variations in fio and h affect al- 
most every aspect of the model, such as the disk mass fraction, 
galaxy ages and circular velocities, star-formation rates, and the 
power-spectrum shape. In Figures [7] and [§] we plot S distributions 
with h = 0.5 and h — 0.75, respectively. Little change is seen 
in the EdS cases, while for low £lo the net result appears to be 
that lower (higher) values of h produce galaxies which are brighter 
(fainter). Changes in h do not, however, strongly affect the scatter. 
Note that the h — 0.5 ACDM model provides reasonable agree- 
ment in K and J, but is not as successful in B or R. The depen- 
dence of the models on Slo can be gleaned from comparing the 
left and right panels of the various Figures. In general, for a fixed 
halo mass, lower values of Slo will yield a larger disk mass and thus 
brighter luminosity. We have seen repeatedly that high values of fio 
invariably predict galaxies which are excessively faint as compared 
to observations. 

In order to gain broader insight into the dependence of our 
results on the cosmological parameters, we generated data for the 
fio = 1 case with V = 0.2 and h = 0.45-0.75, and for fi = 0.2- 
0.5 cases with Y = Q h and h = 0.55-0.75, both for the S dis- 
tribution and for LC distributions with / = 0.5 and 0.75. All were 
COBE normalized. Only a handful of these models produced val- 
ues of p > 0.001 in any band, and those which did invariably had 
values of £loh ~ 0.2, as required by current measurements. EdS 
models generally produce galaxies which are too faint, unless un- 
acceptably low values of h are adopted, and produce too large a 
scatter, unless unacceptably low values of as are adopted. The best 
fitting models, both in terms of p and scatter, were LC distribu- 
tions with COBE normalization, fio ~ 0.3, and h ~ 0.65, sim- 
ilar to the best-fitting values obtained in earlier work (EL96; van 
den Bosch 2000) and in agreement with current estimates. As men- 
tioned above, the S distribution predictions are generally too faint, 
due to the significant fraction of older galaxies arising from the 
high-redshift tail. 



3.3 Spin-Parameter Distribution 

We now turn to the issue of the spin-parameter distribution. Figures 
^| and [jfj correspond to the S and LC distributions from Figures ^ 
and |§J but with the joint probability distribution in A and v now 
taken into account. In the S distributions, this produces virtually no 
effect in K, as expected, but does broaden the distribution of points 
at each mass in the bluer bands. In every case, however the TF scat- 
ter, is roughly similar to, or less than, that in Figure En the broad- 
ening has aligned along the TF relation itself. This arises from the 
fact that for a given mass, higher peaks (which collapse earlier and 



have higher V c ) will have smaller spins leading to higher gas sur- 
face density and thus higher luminosities. A similar result is found 
in the LC case, though here the Zf distribution is weighted toward 
more recent epochs, exaggerating this effect and leading to a larger 
systematic reduction in the scatter, once our color selection crite- 
rion is enforced. Thus, in practice, accounting for the joint prob- 
ability distribution in A and v can reduce the TF scatter by about 
0.15 mag in B to 0.05 mag in K, for plausible cosmogonies. Oth- 
erwise, the spin distribution does not have a great impact on the TF 
relation. As a corollary to this, we find that high-surface-brightness 
and low-surface-brightness disks, as distinguished by different A or 
Rd in our model (see Figure |l4|), are predicted to lie on the same 
TF relation, as is in fact observed (Tully et al. 1998). 

3.4 Chemical Evolution and the IMF 

The role of chemical evolution is investigated in Figure [ll] for the 
flo = 1 S distribution and the ACDM LC distribution. Compared 
with their constant, solar-metallicity counterparts in Figures H and 
these models have TF scatters which are systematically smaller 
by about 0.1 mag in all wavebands. This simply results from the 
fact that relative to constant metallicity systems, older disks (with 
higher V c ) will continually be forming populations of higher metal- 
licity, and thus have a higher integrated luminosity. Conversely, 
younger systems will be relatively dimmer, so that models with an 
evolving metallicity content will naturally align more tightly along 
a luminosity-circular velocity relationship than those with constant 
Z. 

In Figure |l^ we examine the impact of changing the assumed 
IMF. The left panels depict a Salpeter IMF with a = 0.95 (de- 
noted as Salpeter-2), while the right panels represent a Scalo IMF 
(Scalo, 1986), each for the ACDM LC distribution. These IMFs 
both produce relatively more high-mass stars, thus yielding predic- 
tions which are excessively bright in the bluer bands. Note that the 
Salpeter-2 model matches reasonably in K, but not in any other 
band. This is yet another example of the importance of using multi- 
wavelength constraints to discriminate between models which may 
appear to fit in a single band. 



4 A MODEL THAT WORKS 

4.1 The TF Relation at z=0 and z=l 

As an example of what can be accomplished with the theoreti- 
cal framework presented here, we illustrate specifically the results 
for a COBE-normalized ACDM LC distribution with / = 0.5, 
tto ~ 0.3, h = 0.68, a Salpeter IMF, and with the joint dis- 
tribution in A and v, as well as chemical evolution. We refer to 
this particular model as the 'A' model. While this is not necessar- 
ily the best-fitting model over our entire parameter space, it does 
yield remarkable agreement with the data, as shown in Figure [13], 
and these model parameters fall nicely in line with current observa- 
tional constraints. In addition to a and p, the plots list the predicted 
normalizations and slopes [given by a and b in equation (jjj)]. The 
A model appears to fit the slope, normalization, and scatter of the 
TF relation reasonably well in all bands. The predicted scatter in B 
remains about 0.2 mag higher than observed, but there are several 
reasons one might expect this, given our model assumptions. In the 
next section, we show that features such as strong supernova feed- 
back or a flatter density profile towards the core are not necessary 
to successfully predict many of the global spiral-galaxy properties 
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Figure 13. Agreement of the A model. Present-day TF predictions for the A model in B, R, I, and K for a COBE-normalized / = 0.5 LC distribution with 
Qq = 0.3, h = 0.68, a Salpeter IMF, including the joint distribution in A and v as well as chemical evolution. This model, denoted as the A' model, produces 
excellent agreement in all bands (p > 0.10) with roughly the correct amount of scatter (cr ~ 0.4-0.5). 



investigated here. Thus, these features would tend only to produce 
small changes in the amount of gas available for star formation, so 
that their effects would be largely confined to B. Presumably, how- 
ever, their inclusion would result in a slight dimming and tightening 
the B-band predictions, bringing them more in line with the obser- 
vations. 

Figure [w] shows the predictions of the A model for the TF 
relation at z = 1, together with 6 data points from Vogt et al. (1997) 
for spirals in the range 0.5 < z < 1.0. These data have been 
corrected for dust. At z = 1 our model predicts larger values of 
a (i.e., lower zero-points) and steeper slopes for the TF relation in 
every band. This results in the TF relation becoming about 1 mag 



brighter at 50 km s _1 , and about 2 mag brighter at 500 km s _1 
in the B-band, but conspires to produce very little change at these 
scales in the K band. Moreover, the scatter in the TF relation in 
all bands at z = 1 is predicted to be roughly the same as that at 
z = 0, though we do not impose any color-selection criterion in 
the high-redshift case, since it is not as clear what the observational 
selections are. Though the B-band data shown are comprised of 
only 6 points, they have an intrinsic scatter consistent with that of 
the low-z data, and are fit by our model extremely well (p — 0.32), 
though, as in the z — case, the predicted scatter is about 0.2 mag 
too large. 
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Figure 14. High-redshift predictions. TF prediction of the A model for 2 = 1, along with 6 B-band data point from Vogt et al. (1996) for spirals in the range 
0.5 < z < 1. The 2 = 1 prediction have steeper slopes (given by b) and lower zero-points (given by a), than the counterpart z = predictions, resulting in a 
brighter B-band TF relation, but roughly a roughly similar relation in K. 



4.2 The Luminosity Function and Surface-Brightness 
Distribution 

The LF is given by dn/dM x = (dn/dM)dM/dM x where 
dn/dM is the number density of objects, per unit mass, obtained 
from equation (^). An equation for dM/dM\ is derived from the 
model outputs assuming the mass-magnitude relationship obeys a 
power law, which is an excellent assumption at these scales. For 
B-band magnitudes, we apply the dust correction formula of Wang 
& Heckman (1996) in calculating the LF. 

The solid lines in Figure |lj are the predictions for the present- 
day B and K LFs for the A model. The horizontal error bars indi- 
cate the 1-ct error in our mass-magnitude power-law fits. No color 



cuts have been imposed. The dashed line shows the B-band predic- 
tion without correcting for dust. The dot-dashed line is the B-band 
fit to the 2dF Survey data (Folkes et al. 1999), and the open squares 
are data from Zucca et al. (1997). The K-haxiA data are taken from 
Gardner et al. (1997; open triangles) and Glazebrook et al. (1995; 
filled triangles). 

The LF fits are where this model is weakest, as they rely heav- 
ily on the detailed number distribution of galaxies as predicted by 
the PS formalism, which is known to be inaccurate. We find that 
while the model reasonably approximates the faint end of the K- 
band data, it slightly underpredicts the number of galaxies around 
the knee (near 10 Mq) and overpredicts the bright end. We note 
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Figure 15. The luminosity function. B and K LF predictions of the A 
model (solid lines) together with 1-cr errors (horizontal error bars). The 
dashed line is the 73-band prediction before correcting for dust. The dot- 
dashed line is the B-band fit to the 2dF Survey data (Folkes et al. 1999), 
and the open squares are data from Zucca et al. (1997). The iC-band data 
are taken from Gardner et al. (1997; open triangles) and Glazebrook et 
al (1995; filled triangles). 




Figure 16. The surface-brightness-magnitude relation. A-model predic- 
tions for the mean B-band surface brightness as a function of magnitude 
(solid line) at z = 0.4, which corresponds to the mean redshift of the HDF 
data points (open triangles; taken from Driver et al. 1999). The dots are the 
data from Driver & Cross (2000), derived from the 2dF survey. The dashed 
lines show the 1-cr error in the prediction. We have also included the loca- 
tion of the anomalous LSB giant Malin 1 . 



that the bright end of the observed if -band LF is not presently well 
constrained, due to sample size and evolutionary effects. The dust- 
corrected B-band predictions yield reasonable agreement from the 
largest galaxies down to about 10 11 Mq, after which they overshoot 
the faint end by a considerable margin, as PS models are known to 
do. Given the shortcomings and limitations of the standard PS for- 
malism, it is quite remarkable that the model yields as good agree- 
ment as it does over a factor of 100 in both number and luminosity. 



Somerville & Primack (1999), using a modified PS prescription, do 
achieve slightly better agreement. 

Generically, two explanations might resolve the discrepancy at 
the faint end in B, which is seen to persist even when more sophisti- 
cated formalisms than PS theory are implemented. One possibility 
is that our models overestimate the luminosity for a given mass. 
This would be the case, for example, if supernova feedback played 
a significant role in removing the gas. It is unclear to what extent 
the energy they release alters the global properties of spirals. If they 
serve only to "puff up" the gas distribution (Mac Low & Ferrara 
1999), then our model would remain effectively unchanged. If they 
expel gas from the halo, this could alter our predictions. Inspection 
of Figure [l^, however, reveals that, at low masses, feedback would 
have to reduce the predictions by many magnitudes in order to yield 
agreement. Somerville & Primack do find that the inclusion of feed- 
back improves the agreement at the faint end in B, but still returns 
LF values on the high side. Our TF predictions, however, match 
B-band data down to V c = 50 km s _1 . Strong feedback would not 
only disrupt this agreement for low V c , but also alter the predic- 
tions at other wavelengths. Some recent observations do indicate a 
break in the near-IR TF relation at V c ~ 90 km s _1 (McGaugh 
et al. 2000), but these galaxies are found to be extremely gas rich, 
such that feedback is not sufficient to expel the gas. Fine tuning by 
supernovae of the gas mass available for star formation would be 
allowed by our model. This could improve both the TF predictions 
in B and slightly alleviate the LF discrepancy, but it appears un- 
likely that the latter problem can be solved entirely via feedback. 
We point out that the Schmidt law implicitly assumes supernova 
feedback which can redistribute gas within the disk and serves to 
delay star formation. It is only inflows and outflows which are not 
accounted for in our model, and detailed multi-phase ISM simu- 
lations seem to indicate that these are negligible and that spirals 
may be well-approximated as "island universes." Note also that we 
may be implicitly assuming some influence of supernova (or other) 
feedback in keeping the specific angular momentum of the baryons 
equal to that of the halo. 

Another possibility is that the model is simply predicting too 
many small halos. This may be due in part to the breakdown of 
the PS formalism in accurately describing the distribution of ha- 
los over these scales, but may also indicate that the assumed power 
spectrum itself is supplying too much power to these scales. Re- 
cent evidence suggests that current numerical models do indeed 
overpredict the number of dwarf halos surrounding a galaxy such 
as the Milky Way by an order of magnitude or more (Moore et 
al. 2000). Kamionkowski & Liddle (1999) show that this stark dis- 
crepancy can be addressed through plausible modifications to the 
power spectrum. This solution would reconcile the LF predictions 
with observations, while leaving the TF predictions untouched. 

A very interesting test that is usually overlooked is that of 
the surface-brightness-magnitude relation. The solid line in Fig- 
ure |^ shows our A-model prediction for the mean B-band sur- 
face brightness within the disk scale length, Rd [see equation (g)], 
as a function of magnitude; the dashed lines are the 1-a bounds. 
The open triangles are data from Driver et al. (1999), who used a 
volume-limited sample derived from the Hubble Deep Field, and 
the dots are the data from Driver & Cross (2000) derived from the 
2dF survey. We have also included the location of the anomalous 
LSB giant Malin 1. Although we are plotting model predictions for 
z = 0.4, the mean redshift for the galaxies in the HDF data set, the 
brightening of galaxies between z = 0.1 and 0.4 will have a mi- 
nor impact on the outcome of the fit. No color cuts have been im- 
posed. The first thing to note is that the data correspond to effective 
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surface brightness as defined in Driver et al. (1999), and therefore 
include a bulge component which is absent in our computed cen- 
tral surface brightness for an exponential diskf^|. The first striking 
feature of the data itself is that there exists a correlation between 
absolute magnitude and surface brightness, as first discovered by 
Binggeli & Cameron (1991) for the Virgo population and later con- 
firmed for the field population by Driver et al. (1999). Gratifyingly 
enough, our model predicts a similar correlation which produces a 
reasonable fit (albeit slightly fainter for the reasons detailed above 
and with a scatter only 15% larger than that currently observed). 
Until more complete surveys can determine the precise scatter of 
the surface-brightness-magnitude relation, it is too early to con- 
clude if our model predicts a scatter that is in fact too large. As 
already proposed by Driver & Cross (2000), this relation offers a 
physical alternative to the "botanical" Hubble-galaxy classification 
scheme [see Fig. 4 in Driver & Cross (2000)], with the x-axis de- 
termined by the mass and the y-axis determined primarily by the 
angular momentum of the dark halo. Here we confirm that the x- 
axis is indeed determined by the mass, but the y-axis is determined 
both by the spin of the halo (A) and by the redshift of formation 
(i.e., by the initial conditions). It is also reassuring to find that our 
model predicts that LSB giants like Malin- 1 should be very rare ob- 
jects (about 4-<j fluctuations), which agrees well with the observed 
paucity of these objects. As was the case for the TF predictions, 
most other models investigated (with different cosmogonies and so 
on) grossly failed to fit the data. Unlike the TF predictions, how- 
ever, the surface-brightness predictions do depend sensitively on 
A (which fixes the radial scale length), and thus provide a strong 
complementary constraint on the models. 

5 CONCLUSIONS 

We have constructed a largely analytic disk-galaxy model, built 
on that of HJ99, which incorporates such features as the halo 
formation-redshift distribution, the joint distribution in spin and 
peak height, the dependence on the cosmological model, and star 
formation with chemical evolution. Exploring the model's predic- 
tions over a range of wavebands, we find that both parent-halo fea- 
tures and stellar evolution play a role in defining the TF relation 
and other disk-galaxy properties; the former fixes the overall rela- 
tion and spread, while the latter fine tunes these to observed levels. 
Our investigation has yielded a number of specific results: 

• We confirm that the spread in formation redshifts is the pri- 
mary source of scatter in the TF relation. In agreement with pre- 
vious studies, we predict a TF relation which generally broadens 
with decreasing V c and exhibits increasingly large scatter towards 
bluer wavebands (HJ99; van den Bosch 2000; Somerville & Pri- 
mack 1999). 

• We find that Qq — 1 models are generally too faint to match 
the observed TF relation, but contrary to the claim of van den 
Bosch (2000), this result does not depend on the inclusion of adi- 
abatic disk contraction, universal-halo profiles (Navarro, Frenk & 
White 1997), or the precise definition of V c . We find that success- 
ful models tend to favor low values of Q,oh ~ 0.2, and in particular 
f2o ~ 0.3 and h ~ 0.65 (assuming COBE normalization) ,p] 

** Note that the data points are slightly higher than the Freeman law (Free- 
man 1970; Driver et al. 1999). 

tt Recent preliminary results from the BOOMERANG experiment (Lange 
et al. 2000) point to larger values of f2;,/i 2 . A higher baryon fraction would 



• The inclusion of the spin-parameter distribution, with an anti- 
correlation with peak height, tends to act along the TF axis itself 
and can therefore slightly reduce the TF scatter in some models, 
but does not otherwise have a major impact on the TF relation. 
Thus, uncertainties in obtaining accurate A distributions, such as 
nonlinear effects, proper calculation of binding energies, and so on, 
are likely unimportant. The spin distribution is important, however, 
in determining disk scale lengths and surface brightnesses, which 
can provide independent tests of the model. 

• The incorporation of chemical evolution leads to older 
(younger) populations which are brighter (fainter), leading to a re- 
duction of TF scatter of about 0.1 mag in all wavebands. 

• Many models can be found to yield reasonable agreement with 
the data in one or two bands, but not in others. Typically, I and K 
data are the easiest to fit, since these wavelengths are largely de- 
tached from the star-formation history. Imposing multi-wavelength 
constraints, however, provides a strong discriminator among mod- 
els. 

• Models which yield the best agreement with TF data across all 
wavelengths have cosmological-parameter values in good agree- 
ment with current estimates, and require that most disks from in 
the range 0.5 < z < 2.0, with little subsequent accretion, as 
suggested by Peebles (1999). In particular, the reasonable agree- 
ment in B seems to rule out very late collapse. The / = 0.5 LC 
distribution appears to represent a suitable choice for the range of 
halo formation redshifts. Note that this agreement is obtained with 
a spread in Zf of Az//(1 + Zf) ~ 0.5, considerably larger than 
that suggested by the results of Eisenstein & Loeb (1996) and van 
den Bosch (2000), demonstrating the impact of regulating effects 
from stellar evolution and metal enrichment. 

• Models fitting the TF relation at z = also appear to obtain 
the correct, steeper B-band TF relation at z = 1, and predict a 
near-IR z = 1 TF relation which is essentially identical to the local 
one at the scales investigated. 

• Successful models also roughly match the observed surface- 
brightness distribution of spirals and reasonably match the ob- 
served B- and J\~-band LFs. 

• In general, this success argues that there is probably little 
room for stochasticity arising from other complicating effects, such 
as mergers. Spiral galaxies appear to behave like island universes 
(Mac Low & Ferrara 1999; Peebles 1999). Small variations to the 
disk baryon content from such mechanisms could still be tolerated, 
though their effects would be largely confined to the B band. Mod- 
eling these weaker effects might serve to improve the model's pre- 
diction for the B-band TF relation and LF. 

There are several features common to many SAMs which have 
not been incorporated into our framework, as these do not appear 
to strongly impact the properties investigated here. For example, 
we adopt an isothermal profile, rather than the universal form of 
Navarro, Frenk, & White (1997). While observations, as well as 
some simulations, indicate that real halos do resemble truncated 
isothermal spheres, rather than universal profiles (Sellwood 1999; 
Dubinski & Carlberg 1991; Spergel & Hernquist 1992), our model 
is clearly too steep in the core. For our purposes, however, changes 
to the profile of a halo with a fixed mass are only relevant inas- 

brighten the TF predictions, but this could be counteracted by a slightly 
higher value of h (or by changing other input parameters so as to still 
remain within current constraints), so as to yield a similarly well-fitting 
model. Many of the conclusions presented here are robust with respect to 
the precise values of these various parameters. 



© 0000 RAS, MNRAS 000, 000-000 



Galacto synthesis: Halo Histories, Star Formation, and Disks 23 



much as they alter the local gas surface density and thus the SFR. 
This is similar to the effect of changing the spin parameter, and we 
have seen that our model changes relatively little over the broad 
range investigated in A and thus E (HJ99). Jimenez et al. (1997) 
developed a similar disk-galaxy model, but explicitly incorporating 
the universal profile, and found little difference from the results of 
HJ99. The behavior of an isothermal profile at the core, however, 
will have the effect of making younger galaxies too bright in B. As 
discussed above, this might partially account for the slightly poorer 
TF fit and larger predicted scatter in B, as well as the discrepancy 
with the S-band LF. 

While the assumption of monolithic disk collapse cannot be 
strictly accurate, accounting for adiabatic disk contraction also does 
not appear essential to predicting these disk properties. Moreover, 
adiabatic contraction tends to make spirals more dark-matter dom- 
inated at the core, whereas observations indicate that they are more 
baryon dominated (Sellwood 1999). Our model depends on an ac- 
curate assessment of V c , and not on the details of disk formation, 
which do not affect our ultimate results. Similarly, the inclusion 
of a stability threshold criterion for star formation does not appear 
essential to our results, and some authors have argued against the 
importance of such a threshold (Ferguson et al. 1998). 

Another key assumption in this model is that the specific angu- 
lar momentum of the baryons is the same as that of the dark matter. 
This is a common assumption of analytic calculations (e.g. Fall & 
Efstathiou 1980; Mo, Mao & White 1998), but it has been argued 
from hydrodynamic simulations that this is not the case; cooling of 
gas leads to substructure which can couple with the halo, leading 
to loss of a sizeable fraction of the baryon angular momentum (e.g. 
Navarro & Benz 1991; Weil, Eke & Efstathiou; Navarro & Stein- 
metz 2000). The predictions of our model would thus be altered 
significantly, with much smaller disks, higher early star formation 
rates and different colors. Therefore we need to assume that such 
angular momentum loss does not occur in practice. It is fair to say 
that the question of how (or indeed if) baryon angular momentum 
is conserved is open at present, although there are plenty of ideas 
for how the coupling can be prevented. These include gas ejection 
from sub-units by supernovae (Efstathiou 2000), suppression of gas 
cooling until the halo is established with a smooth profile (Weil, 
Eke & Efstathiou 1998), prevention of formation of small halos by 
a cut-off in the power spectrum (Kamionkowski & Liddle 1999), or 
changing the nature of the dark matter (Spergel & Steinhardt 1999). 
In addition, van Kampen (2000) proposes that numerical simula- 
tions underestimate dynamical friction of the sub-units, leading to 
overestimated angular-momentum transfer. 

A generic problem is the possibility that our model either 
over-predicts or under-predicts the amount of gas which is con- 
verted into disk stars, due to effects such as supernova feedback, 
the presence of hot x-ray gas, or enhanced star formation from late 
infall. Though we have assumed essentially featureless (in space 
and time) star formation, the addition or removal of large amounts 
of gas via such stochastic mechanisms would increase the scatter in 
the TF relation and shatter the agreement obtained here. To the ex- 
tent that such processes result in smaller fluctuations to the amount 
of gas available for star formation, their effects would be confined 
primarily to B, where our models are indeed weakest. Mergers 
likely play only a very limited role. It is well known that classic 
spirals could not maintain their thin-disk structure in the case of 
extreme merging. EL96 found that excluding objects which had 
merged with a mass > 20% (as opposed to 50%) of the parent 
halo, only reduced the predicted J-band scatter by 20%, suggest- 
ing that accretion does not play a major role, or possibly induces 



scatter along the near-infrared TF relation. Either way, this would 
imply that the TF relation at higher redshift looks similar to that 
at z = 0, and this is precisely what we find when comparing the 
z — and z = 1 TF relations in / and K. Future work will be 
aimed at quantitatively assessing the precise impact of these vari- 
ous other mechanisms, as well as obtaining predictions for the high- 
redshift universe. While our model is a long way from providing a 
complete picture of galactosynthesis, our results suggest that future 
models addressing spiral properties must invariably incorporate de- 
tailed star formation, as well as initial conditions, and that features 
such as exponential disks, little or weak substructure, Schmidt-law 
star formation, and little gas inflow or outflow, should be generic 
predictions of more sophisticated modelling. 
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APPENDIX 

The luminosity of an arbitrary stellar population can be computed 
analytically as follows. Simple stellar populations (SSPs) are the 
building blocks of any arbitrarily complicated population since the 
latter can be computed as a sum of SSPs, once the star-formation 
rate is provided. In other words, the luminosity of a stellar popula- 
tion of age to (since the beginning of star formation), in waveband 
A, can be written as: 
/•to rZf 

L x (t )= / / L SS p,\(Z,to-t)dZdt (11) 

J0 J Zi 

where the luminosity of the SSP is: 

Lssp,\{Z, to — t) = 
rM max 

/ SFR(Z,M,t)l x (Z,M,t -t) — dM (12) 

and l\(Z, M, to —t) is the luminosity of a star of mass M, metallic- 
ity Z, and age to— t, Zi and Zf are the initial and final metallicities, 
dN/dM represents the IMF, M m ;„ and M max are the smallest and 
largest stellar mass in the population, and SFR(Z, M, t) is the star 
formation rate at the time t when the SSP is formed. 

The magnitudes for a SSP (normalized to 1 Mq) as a function 
of age and metallicity, for a given photometric band UBVRIJK, are 
approximated to within 4% by: 

4 4 

M x = -2.5 x Yl X ' Cx(i + + 1)YJ ' (13) 

i=0 j=0 

where 

X = 5.76 + 3.18 log r + 1.26 log 2 r + 2.64 log 3 r 

+1.81 log 4 r + 0.38 log 5 r, (14) 

Y = 2.0 + 2.059 log C + 1.041 log 2 < 

+0.172 log 3 C- 0.042 log 4 C, (15) 
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and 



t 

Gyr' 



«-f- 



(16) 



Luminosities 

are obtained simply from L A = iO" - 4 *( M ©A-Af A ) 5 where 
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nents of X and Y, respectively, and as indices defining elements of 
the C\ matrices, given by 
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